Safety factor diagnostic for tokamak core plasma from three-dimensional reconstruction of pellet ablation trail.

Using a stereo camera system, a new diagnostic for the safety factor of the core plasma based on the pellet ablation trail is applied on the Experimental Advanced Superconducting Tokamak (EAST). In EAST discharge No. 128 874, a shattered pellet injection system is applied to inject a shattered neon pellet into the EAST. Since the strong magnetic field in tokamaks binds the ablated pellet material, the orientation of the pellet ablation trail is the same as the local magnetic field direction. Thus, from the three-dimensional reconstruction result of the pellet ablation trail, the local safety factor q can be obtained. The motional Stark effect (MSE) diagnostic is applied to determine the safety factor q profile in this shot. The determined safety factor q results for this new diagnostic are in quantitative agreement with those from the MSE diagnostic with the mean relative difference of only 6.8%, confirming the effectiveness of this new diagnostic of the safety factor.

Life expectancy, long-term care demand and dynamic financing mechanism simulation: an empirical study of Zhejiang Pilot, China.

BACKGROUND: China has piloted Long-Term Care Insurance (LTCI) to address increasing care demand. However, many cities neglected adjusting LTCI premiums since the pilot, risking the long-term sustainability of LTCI. Therefore, using Zhejiang Province as a case, this study simulated mortality-adjusted long-term care demand and the balance of LTCI funds through dynamic financing mechanism under diverse life expectancy and disability scenarios. METHODS: Three-parameter log-quadratic model was used to estimate the mortality from 1990 to 2020. Mortality with predicted interval from 2020 to 2080 was projected by Lee-Carter method extended with rotation. Cohort-component projection model was used to simulate the number of older population with different degrees of disability. Disability data of the older people is sourced from China Health and Retirement Longitudinal Study 2018. The balance of LTCI fund was simulated by dynamic financing actuarial model. RESULTS: Life expectancy of Zhejiang for male (female) is from 80.46 (84.66) years in 2020 to 89.39 [86.61, 91.74] (91.24 [88.90, 93.25]) years in 2080. The number of long-term care demand with severe disability in Zhejiang demonstrates an increasing trend from 285 [276, 295] thousand in 2023 to 1027 [634, 1657] thousand in 2080 under predicted mean of life expectancy. LTCI fund in Zhejiang will become accumulated surplus from 2024 to 2080 when annual premium growth rate is 5.25% [4.20%, 6.25%] under various disability scenarios, which is much higher than the annual growth of unit cost of long-term care services (2.25%). The accumulated balance of LTCI fund is sensitive with life expectancy. CONCLUSIONS: Dynamic growth of LTCI premium is essential in dealing with current deficit around 2050 and realizing Zhejiang's LTCI sustainability in the long-run. The importance of dynamic monitoring disability and mortality information is emphasized to respond immediately to the increase of premiums. LTCI should strike a balance between expanding coverage and controlling financing scale. This study provides implications for developing countries to establish or pilot LTCI schemes.

CHANGE-seq-BE enables simultaneously sensitive and unbiased in vitro profiling of base editor genome-wide activity.

Base editors ( BE ) enable programmable conversion of nucleotides in genomic DNA without double-stranded breaks and have substantial promise to become new transformative genome editing medicines. Sensitive and unbiased detection of base editor off-target effects is important for identifying safety risks unique to base editors and translation to human therapeutics, as well as accurate use in life sciences research. However, current methods for understanding the global activities of base editors have limitations in terms of sensitivity or bias. Here we present CHANGE-seq-BE, a novel method to directly assess the off-target profile of base editors that is simultaneously sensitive and unbiased. CHANGE-seq-BE is based on the principle of selective sequencing of adenine base editor modified genomic DNA in vitro , and provides an accessible, rapid, and comprehensive method for identifying genome-wide off-target mutations of base editors.

Integration of Multi-Omics, Histological, and Biochemical Analysis Reveals the Toxic Responses of Nile Tilapia Liver to Chronic Microcystin-LR Exposure.

Microcystin-LR (MC-LR) is a cyanobacterial metabolite produced during cyanobacterial blooms and is toxic to aquatic animals, and the liver is the main targeted organ of MC-LR. To comprehensively understand the toxicity mechanism of chronic exposure to environmental levels of MC-LR on the liver of fish, juvenile Nile tilapia were exposed to 0 µg/L (control), 1 µg/L (M1), 3 µg/L (M3), 10 µg/L (M10), and 30 µg/L (M30) MC-LR for 60 days. Then, the liver hepatotoxicity induced by MC-LR exposure was systematically evaluated via histological and biochemical determinations, and the underlying mechanisms were explored through combining analysis of biochemical parameters, multi-omics (transcriptome and metabolome), and gene expression. The results exhibited that chronic MC-LR exposure caused slight liver minor structural damage and lipid accumulation in the M10 group, while resulting in serious histological damage and lipid accumulation in the M30 group, indicating obvious hepatotoxicity, which was confirmed by increased toxicity indexes (i.e., AST, ALT, and AKP). Transcriptomic and metabolomic analysis revealed that chronic MC-LR exposure induced extensive changes in gene expression and metabolites in six typical pathways, including oxidative stress, apoptosis, autophagy, amino acid metabolism, primary bile acid biosynthesis, and lipid metabolism. Taken together, chronic MC-LR exposure induced oxidative stress, apoptosis, and autophagy, inhibited primary bile acid biosynthesis, and caused fatty deposition in the liver of Nile tilapia.

Global multifaceted biodiversity patterns, centers, and conservation needs in angiosperms.

The Convention on Biological Diversity seeks to conserve at least 30% of global land and water areas by 2030, which is a challenge but also an opportunity to better preserve biodiversity, including flowering plants (angiosperms). Herein, we compiled a large database on distributions of over 300,000 angiosperm species and the key functional traits of 67,024 species. Using this database, we constructed biodiversity-environment models to predict global patterns of taxonomic, phylogenetic, and functional diversity in terrestrial angiosperms and provide a comprehensive mapping of the three diversity facets. We further evaluated the current protection status of the biodiversity centers of these diversity facets. Our results showed that geographical patterns of the three facets of plant diversity exhibited substantial spatial mismatches and nonoverlapping conservation priorities. Idiosyncratic centers of functional diversity, particularly of herbaceous species, were primarily distributed in temperate regions and under weaker protection compared with other biodiversity centers of taxonomic and phylogenetic facets. Our global assessment of multifaceted biodiversity patterns and centers highlights the insufficiency and unbalanced conservation among the three diversity facets and the two growth forms (woody vs. herbaceous), thus providing directions for guiding the future conservation of global plant diversity.

Interspecies regulatory landscapes and elements revealed by novel joint systematic integration of human and mouse blood cell epigenomes.

Knowledge of locations and activities of cis -regulatory elements (CREs) is needed to decipher basic mechanisms of gene regulation and to understand the impact of genetic variants on complex traits. Previous studies identified candidate CREs (cCREs) using epigenetic features in one species, making comparisons difficult between species. In contrast, we conducted an interspecies study defining epigenetic states and identifying cCREs in blood cell types to generate regulatory maps that are comparable between species, using integrative modeling of eight epigenetic features jointly in human and mouse in our V al i dated S ystematic I ntegrati on (VISION) Project. The resulting catalogs of cCREs are useful resources for further studies of gene regulation in blood cells, indicated by high overlap with known functional elements and strong enrichment for human genetic variants associated with blood cell phenotypes. The contribution of each epigenetic state in cCREs to gene regulation, inferred from a multivariate regression, was used to estimate epigenetic state Regulatory Potential (esRP) scores for each cCRE in each cell type, which were used to categorize dynamic changes in cCREs. Groups of cCREs displaying similar patterns of regulatory activity in human and mouse cell types, obtained by joint clustering on esRP scores, harbored distinctive transcription factor binding motifs that were similar between species. An interspecies comparison of cCREs revealed both conserved and species-specific patterns of epigenetic evolution. Finally, we showed that comparisons of the epigenetic landscape between species can reveal elements with similar roles in regulation, even in the absence of genomic sequence alignment.

Comprehensive mRNA and microRNA analysis revealed the effect and response strategy of freshwater fish, grass carp (Ctenopharyngodon idella) under geosmin exposure.

Geosmin is an environmental pollutant that causes off-flavor in water and aquatic products. The high occurrence of geosmin contamination in aquatic systems and aquaculture raises public awareness, however, few studies have investigated the response pathways of geosmin stress on freshwater fish. In this research, grass carp were exposed to 50 µg/L geosmin for 96 h, liver tissue was sequenced and validated using real-time qPCR. In total of 528 up-regulated genes and 488 down-regulated genes were observed, includes cytochrome P450 and uridine diphosphate (UDP)-glucuronosyltransferase related genes. KEGG analysis showed that chemical carcinogenesis-DNA adducts, metabolism of xenobiotics by cytochrome P450, drug metabolism-cytochrome P450 pathway was enriched. Common genes from the target genes of microRNAs and differential expression genes are enriched in metabolism of xenobiotics cytochrome P450 pathway. Two miRNAs (dre-miR-146a and miR-212-3p) down regulated their target genes (LOC127510138 and adh5, respectively) which are enriched cytochrome P450 related pathway. The results present that geosmin is genetoxic to grass carp and indicate that cytochrome P450 system and UDP-glucuronosyltransferase play essential roles in biotransformation of geosmin. MicroRNAs regulate the biotransformation of geosmin by targeting specific genes, which contributes to the development of strategies to manage its negative impacts in both natural and artificial environments.

Comprehensive Gas Prevention and Control Technique for Mining the First Seam in Short-Distance Outburst Coal Seam Groups.

Coal and gas outbursts are a phenomenon whereby broken coal and gas suddenly erupt from the coal body into the mining space under pressure. The Diandong mining area is a group of close-range outburst coal seams in which the gas content is up to 20 m3/t and gas pressure can reach 3 MPa. Research has been conducted on engineering challenges such as advanced detection and prevention of interlayer excavation in close-range coal seam groups, improvement of gas extraction quality and efficiency in low-permeability coal seam groups, and traceability and evaluation of joint extraction of coal seam groups. Through this study, advanced detection technology with full coverage in front of the excavation working face has been constructed as well as advanced pre-extraction technology for adjacent coal seams and this coal seam in ultraclose layers. We have developed a method for achieving the standard of cross-layer fixed-point hole expansion and permeability enhancement for the first mining of coal seams in a coal seam group. A combined process of graded enhanced pre-extraction and segmented regulation and extraction was proposed, which included "fixed-point control section sealing pre-extraction of coal seam groups and secondary sealing and extraction of mining pressure relief orifice."

Simulation and Experimental of Infiltration and Solidification Process for Al2O3(3D)/5083Al Interpenetrating Phase Composite for High Speed Train Prepared by Low-Pressure Infiltration.

Understanding the infiltration and solidification processes of liquid 5083Al alloy into Al2O3 three-dimensional reticulated porous ceramic (Al2O3(3D) RPC) is essential for optimizing the microstructure and properties of Al2O3(3D)/5083Al interpenetrating phase composites (IPCs) prepared by low-pressure infiltration process (LPIP). This study employs ProCAST software to simulate the infiltration and solidification processes of liquid 5083Al with pouring velocities (PV) of 0.4 m/s infiltrating into Al2O3(3D) RPC preforms with varying porosities at different pouring temperatures (PT) to prepare Al2O3(3D)/5083Al IPCs using LPIP. The results demonstrate that pore diameter of Al2O3(3D) RPC preforms and PT of liquid 5083Al significantly influence the of the infiltration. Solidification process analysis reveals that the Al2O3(3D) RPC preform with smaller pore diameters allows the lower pouring velocity of 5083Al to solidify faster compared to the preform with larger pore diameters. Al2O3(3D)/5083Al IPCs were prepared successfully from Al2O3(3D) RPC porosity of 15 PPI with liquid 5083Al at PV 0.4 m/s and PT 800 °C using LPIP, resulting in nearly fully dense composites, where both Al2O3(3D) RPCs and 5083Al interpenetrate throughout the microstructure. The infiltration and solidification defects were reduced under air pressure of 0.3 MPa (corresponding to PV of 0.4 m/s) during LPIP. Finite volume method simulations are in good agreement with experimental data, validating the suitability of the simplified model for Al2O3(3D) RPCs in the infiltration simulation.

Global distribution and evolutionary transitions of floral symmetry in angiosperms.

Floral symmetry plays an important role in plant-pollinator interactions and may have remarkable impacts on angiosperm diversification. However, spatiotemporal patterns in floral symmetry and drivers of these patterns remain unknown. Here, using newly compiled floral symmetry (actinomorphy versus zygomorphy) data of 279,877 angiosperm species and their distributions and phylogenies, we estimated global geographic patterns and macroevolutionary dynamics of floral symmetry. We found that frequency of actinomorphic species increased with latitude, while that of zygomorphic species decreased. Solar radiation, present-day temperature, and Quaternary temperature change correlated with geographic variation in floral symmetry frequency. Evolutionary transitions from actinomorphy to zygomorphy dominated floral symmetry evolution, although the transition rate decreased with decreasing paleotemperature throughout the Cenozoic. Notably, we found that zygomorphy may not favor diversification of angiosperms as previously observed in some clades. Our study demonstrates the influence of (paleo)climate on spatiotemporal patterns in floral symmetry and challenges previous views about role of flower symmetry in angiosperm diversification.

Potent and uniform fetal hemoglobin induction via base editing.

Inducing fetal hemoglobin (HbF) in red blood cells can alleviate ß-thalassemia and sickle cell disease. We compared five strategies in CD34+ hematopoietic stem and progenitor cells, using either Cas9 nuclease or adenine base editors. The most potent modification was adenine base editor generation of γ-globin -175A>G. Homozygous -175A>G edited erythroid colonies expressed 81 ± 7% HbF versus 17 ± 11% in unedited controls, whereas HbF levels were lower and more variable for two Cas9 strategies targeting a BCL11A binding motif in the γ-globin promoter or a BCL11A erythroid enhancer. The -175A>G base edit also induced HbF more potently than a Cas9 approach in red blood cells generated after transplantation of CD34+ hematopoietic stem and progenitor cells into mice. Our data suggest a strategy for potent, uniform induction of HbF and provide insights into γ-globin gene regulation. More generally, we demonstrate that diverse indels generated by Cas9 can cause unexpected phenotypic variation that can be circumvented by base editing.

ETV6 represses inflammatory response genes and regulates HSPC function during stress hematopoiesis in mice.

ETS variant 6 (ETV6) encodes a transcriptional repressor expressed in hematopoietic stem and progenitor cells (HSPCs), where it is required for adult hematopoiesis. Heterozygous pathogenic germline ETV6 variants are associated with thrombocytopenia 5 (T5), a poorly understood genetic condition resulting in thrombocytopenia and predisposition to hematologic malignancies. To elucidate how germline ETV6 variants affect HSPCs and contribute to disease, we generated a mouse model harboring an Etv6R355X loss-of-function variant, equivalent to the T5-associated variant ETV6R359X. Under homeostatic conditions, all HSPC subpopulations are present in the bone marrow (BM) of Etv6R355X/+ mice; however, these animals display shifts in the proportions and/or numbers of progenitor subtypes. To examine whether the Etv6R355X/+ mutation affects HSPC function, we performed serial competitive transplantation and observed that Etv6R355X/+ lineage-sca1+cKit+ (LSK) cells exhibit impaired reconstitution, with near complete failure to repopulate irradiated recipients by the tertiary transplant. Mechanistic studies incorporating cleavage under target and release under nuclease assay, assay for transposase accessible chromatin sequencing, and high-throughput chromosome conformation capture identify ETV6 binding at inflammatory gene loci, including multiple genes within the tumor necrosis factor (TNF) signaling pathway in ETV6-sufficient mouse and human HSPCs. Furthermore, single-cell RNA sequencing of BM cells isolated after transplantation reveals upregulation of inflammatory genes in Etv6R355X/+ progenitors when compared to Etv6+/+ counterparts. Corroborating these findings, Etv6R355X/+ HSPCs produce significantly more TNF than Etv6+/+ cells post-transplantation. We conclude that ETV6 is required to repress inflammatory gene expression in HSPCs under conditions of hematopoietic stress, and this mechanism may be critical to sustain HSPC function.

Combined effects of dietary faba bean water extract and vitamin K3 on growth performance, textural quality, intestinal characteristics, oxidative and immune responses in grass carp.

Faba bean water extract (FBW) and vitamin K3 (VK3) have been demonstrated to improve the muscle textural quality of fish. To better apply these two feed additives in commercial aquaculture setting, four experimental diets (control, commercial feed group; 15% FBW, 15% faba bean water extract group; 2.5% VK3, 2.5% vitamin K3 group; combined group, 15% faba bean water extract + 2.5% vitamin K3 group) were formulated to explore their combined effects of FBW and VK3 on the growth, health status, and muscle textural quality of grass carp. The growth performance, textural quality, intestinal characteristics, and oxidative and immune responses were analyzed on days 40, 80 and 120. The results showed that supplementation with higher doses of FBW and VK3 have no influence on growth-related parameters and immune parameters of grass carp. Notably, compared with the control, fish in the combined group had the highest textural qualities (hardness, chewiness and adhesiveness), followed by those in 15% FBW and 2.5% VK3 groups (P < 0.05). Also, FBW and VK3, to some extent, may lower antioxidative ability of grass carp, as illustrated by lower levels of GSH and CAT in 15% FBW, 2.5% VK3, and combined groups on day 120 (P < 0.05). In addition, enhanced lipase activity was observed in the 15% FBW group. Taken together, the combined supplementation of FBW and VK3 was demonstrated to be a more advanced option than their individual supplementation in a commercial setting owing to the resulting combined effects on both the textural quality and health status of grass carp.

Effects of heat stress on the chemical composition, oxidative stability, muscle metabolism, and meat quality of Nile tilapia (Oreochromis niloticus).

The present study investigated the effects of chronic heat stress (HS) on the chemical composition, oxidative stability, muscle metabolism, and meat quality of Nile tilapia (Oreochromis niloticus). Compared with the control (26 °C), chronic HS (32 °C) lowered growth performance, the contents of whole-body lipid, muscle protein, and muscle lipid. Also, HS significantly increased the contents of reactive oxygen species (ROS) and decreased antioxidative status, causing a decline in meat quality, including increased lipid and protein oxidation, the centrifugal water loss, and cooking loss as well as decreased the fragmentation index and pH at 24 h, which may be attributed to induced apoptosis by excessive ROS in Nile tilapia meat. Moreover, metabolomic analysis showed HS lowered flavor and nutritional value by affecting amino acid, lipid, and nucleotide metabolism. These results reveal that HS adversely affects oxidative stability, meat quality, flavor, and nutrition, warranting its recognition and prevention.

Multiplex-GAM: genome-wide identification of chromatin contacts yields insights overlooked by Hi-C.

Technology for measuring 3D genome topology is increasingly important for studying gene regulation, for genome assembly and for mapping of genome rearrangements. Hi-C and other ligation-based methods have become routine but have specific biases. Here, we develop multiplex-GAM, a faster and more affordable version of genome architecture mapping (GAM), a ligation-free technique that maps chromatin contacts genome-wide. We perform a detailed comparison of multiplex-GAM and Hi-C using mouse embryonic stem cells. When examining the strongest contacts detected by either method, we find that only one-third of these are shared. The strongest contacts specifically found in GAM often involve 'active' regions, including many transcribed genes and super-enhancers, whereas in Hi-C they more often contain 'inactive' regions. Our work shows that active genomic regions are involved in extensive complex contacts that are currently underestimated in ligation-based approaches, and highlights the need for orthogonal advances in genome-wide contact mapping technologies.

Ex vivo prime editing of patient haematopoietic stem cells rescues sickle-cell disease phenotypes after engraftment in mice.

Sickle-cell disease (SCD) is caused by an A·T-to-T·A transversion mutation in the ß-globin gene (HBB). Here we show that prime editing can correct the SCD allele (HBBS) to wild type (HBBA) at frequencies of 15%-41% in haematopoietic stem and progenitor cells (HSPCs) from patients with SCD. Seventeen weeks after transplantation into immunodeficient mice, prime-edited SCD HSPCs maintained HBBA levels and displayed engraftment frequencies, haematopoietic differentiation and lineage maturation similar to those of unedited HSPCs from healthy donors. An average of 42% of human erythroblasts and reticulocytes isolated 17 weeks after transplantation of prime-edited HSPCs from four SCD patient donors expressed HBBA, exceeding the levels predicted for therapeutic benefit. HSPC-derived erythrocytes carried less sickle haemoglobin, contained HBBA-derived adult haemoglobin at 28%-43% of normal levels and resisted hypoxia-induced sickling. Minimal off-target editing was detected at over 100 sites nominated experimentally via unbiased genome-wide analysis. Our findings support the feasibility of a one-time prime editing SCD treatment that corrects HBBS to HBBA, does not require any viral or non-viral DNA template and minimizes undesired consequences of DNA double-strand breaks.

Numerical investigations of the micro lunar dust particles deposition in the human oral respiratory airway.

Understanding the deposition of lunar dust (LD) particles in the human respiratory system is of great significance for protecting astronauts' health from the toxicity of lunar dust. A Euler-Lagrangian approach is adopted to track the LD particle motion in a human oral airway model. The investigations are conducted considering different inspiration rates and micro-particle sizes as well as different abnormal pressures and abnormal temperatures. It is found that 1) almost all the LD particles tend to enter the right lung rather than the left lung, especially in the upper right lobe; 2) at lower ambient pressure, fewer LD particles will deposit in the upper airway, while more particles will enter the lung; 3) at lower temperature, more LD particles are deposited in the upper airway, while fewer are deposited in the lung. In summary, the present work has shown that the LD particles have different depositing properties in the upper airway and the lung lobe regions up to the particle size, inspiration flow rate, temperature and pressure. It should pay more attentions on the upper airway and right upper lobe when it studies the toxicity of the lunar dust, and can't ignore the impact of the environmental temperature and pressure.

An NFIX-mediated regulatory network governs the balance of hematopoietic stem and progenitor cells during hematopoiesis.

The transcription factor (TF) nuclear factor I-X (NFIX) is a positive regulator of hematopoietic stem and progenitor cell (HSPC) transplantation. Nfix-deficient HSPCs exhibit a severe loss of repopulating activity, increased apoptosis, and a loss of colony-forming potential. However, the underlying mechanism remains elusive. Here, we performed cellular indexing of transcriptomes and epitopes by high-throughput sequencing (CITE-seq) on Nfix-deficient HSPCs and observed a loss of long-term hematopoietic stem cells and an accumulation of megakaryocyte and myelo-erythroid progenitors. The genome-wide binding profile of NFIX in primitive murine hematopoietic cells revealed its colocalization with other hematopoietic TFs, such as PU.1. We confirmed the physical interaction between NFIX and PU.1 and demonstrated that the 2 TFs co-occupy super-enhancers and regulate genes implicated in cellular respiration and hematopoietic differentiation. In addition, we provide evidence suggesting that the absence of NFIX negatively affects PU.1 binding at some genomic loci. Our data support a model in which NFIX collaborates with PU.1 at super-enhancers to promote the differentiation and homeostatic balance of hematopoietic progenitors.

Two new sesquiterpenoids from plant endophytic fungus Flammulina velutipes.

Two new sesquiterpenoids, flammupin A (1) and flammupin B (2), along with two known compounds, enokipodin C (3) and 5,5'-dibuthoxy-2,2'-bifuran (4) were obtained from Flammulina velutipes, an endophytic fungus isolated from the roots of Caulophyllum robustum Maxim. The structures were elucidated by the combination of HR-ESI-MS, NMR, and ECD analyses. Compound 3 exhibited moderate to potent cytotoxicity against A549, HeLa, and SMMC-7721 cells with IC50 values ranged from 3.69 to 11.84 µM.

Comparative genome-wide methylation analysis reveals epigenetic regulation of muscle development in grass carp (Ctenopharyngodon idellus) fed with whole faba bean.

Crisp grass carp (CGC), the most representative improved varieties of grass carp (GC), features higher muscle hardness after feeding faba bean (Vicia faba L.) for 90-120 days. DNA methylation, a most widely studied epigenetic modification, plays an essential role in muscle development. Previous studies have identified numerous differentially expressed genes (DEGs) between CGC and GC. However, it remains unknown if the expression levels of these DEGs are influenced by DNA methylation. In the present study, we performed a comprehensive analysis of DNA methylation profiles between CGC and GC, and identified important candidate genes related to muscle development coupled with the transcriptome sequencing data. A total of 9,318 differentially methylated genes (DMGs) corresponding to 155,760 differentially methylated regions (DMRs) were identified between the two groups under the CG context in promoter regions. Combined with the transcriptome sequencing data, 14 key genes related to muscle development were identified, eight of which (gsk3b, wnt8a, wnt11, axin2, stat1, stat2, jak2, hsp90) were involved in muscle fiber hyperplasia, six of which (tgf-ß1, col1a1, col1a2, col1a3, col4a1, col18a1) were associated with collagen synthesis in crisp grass carp. The difference of methylation levels in the key genes might lead to the expression difference, further resulting in the increase of muscle hardness in crisp grass carp. Overall, this study can help further understand how faba bean modulates muscle development by the epigenetic modifications, providing novel insights into the texture quality improvement in other aquaculture fish species by nutritional programming.